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g ■ 1 Problem 
o 

CN I Dirac has written "Each photon then interferes only with itself. Interference between 

O ■ two different photons never occurs." Indeed, a practical definition is that "classical" optics 

Q ■ consists of phenomena due to the interference of photons only with themselves. However, 

photons obey Bose statistics which implies a "nonclassical" tendency for them to "bunch" . 

For a simple example of nonclassical optical behavior, consider two pulses containing 
ni and n2 photons of a single frequency that are simultaneously incident on two sides of a 
lossless, 50:50 beam splitter, as shown in the figure. Deduce the probability that A'^i photons 
■ are observed in the direction of beam 1, where < A'^i < ni + ^2 for a lossless splitter. 

O 

(N 

o 

> 

X 

^ reflected photon {i.e., of the reflected wave from a single input beam) is 90° different from 

that of a transmitted photon |2j. Consider first the cases that rii or n2 is zero. 
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Hint: a relatively elementary argument can be given by recalling that the phase of a 



2 Solution 

An elegant solution can be given by noting the the creation and annihilation operators 
relevant to a beam splitter obey an SU(2) symmetry Here, we give a more elementary 

solution, in the spirit of Feynman j^j . 

Experimental demonstration of the case where ni = 77,2 = 1 was first given in jB] , and the 
case of 77-1 = 77.2 = 2 has been studied in [Zj. 
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2.1 A Single Input Beam 

We first consider the case of a single input beam witli ni > 0. Tlien, of course, 77.2 = 0. 

In a classical view, the input beam would have energy ui = riihuj, where u is the angular 
frequency of the photons. Then, the effect of the 50:50 beam splitter would be to create 
output beams of equal energies, Ui = U2 = ui/2. In terms of photon numbers, the classical 
view would imply that the only possibility for the output beams is Ni = N2 = ni/2. 

But in fact, the transmitted beam can contain any number A'^i of photons between and 
rzi, while the reflected beam contains N2 — ni — Ni photons. 

If the photons were distinguishable, we would assign a probability of (1/2)"^ to each 
configuration of transmitted and reflected photons in the 50:50 splitter. But the photons are 
indistinguishable, so that the probability that A'^i out of ni photons are transmitted is larger 
than (1/2)"^ by the number of ways the rii photons can be arranged into a group of Ni 
transmitted and rii — Ni reflected photons with regard to their order, i.e., by the binomial 
coefficient, 

rill 



ni 



C 

Ni\{ni-Ni)\' 

Thus, the probability P{Ni, ni — Ni\ni, 0) that A^^i out of ni photons (in a single input beam) 
are transmitted by the beam splitter is 

P(iVi,ni-iVi|ni,0) = C;^;Q"\ (2) 

The result is already very nonclassical, in that there is a small, but nonzero probability 
that the entire input beam is transmitted, or reflected. However, in the limit of large rii the 
largest probability is that the numbers of photons in the reflected and transmitted beams 
are very nearly equal. We confirm this by use of Stirling's approximation for large n, 



n\ ^ e"VV27m. (3) 
For large n, and k = {1 + e)n/2, we have 
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(4) 



2n+l 



v/2^(l + ne2/2) ' 

The probability of k photons out of n being transmitted drops to 1/2 the peak probability 



when e ~ yS/n. Hence, for large n the number distribution of photons in the transmitted 
(and reflected) beam is essentially a delta function centered at ri/2, in agreement with the 
classical view. 

The most dramatic difference between the classical and quantum behavior of a single 
beam in a 50:50 beam splitter occurs when ni = 2, 

P(0,2|2,0) = ^, P(l,l|2,0) = ^, P(2,0|2,0) = ^. (5) 
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In the subsequent analysis we shall need to consider interference effects, so we note that 
the magnitude of the probability amplitude that k out of n photons in a single beam are 
transmitted by a 50:50 beam splitter can obtained by taking the square root of eq. Q, 

\A{k,n-k\n,0)\ = ^(^^Y\ (6) 

These amplitudes have the obvious symmetries, 

\A{k,n- k\n,0)\ = \A{n - k, k\n,0)\ = \A{k,n - k\n,0)\ = \A{n - k, k\n,0)\ . (7) 

We must also consider the phases of these amplitudes, or at least the relative phases. The 
hint is that we may consider the phase of a reflected photon to be shifted with respect to that 
of a transmitted photon by 90°, as follows from a classical analysis of waves in a 50:50 beam 
splitter [2] (see also the Appendix). In this problem, we define the phase of a transmitted 
photon to be zero, so that the probability amplitude should include a factor of i = a/— 1 for 
each reflected photon. Thus, we have 



A{k,n-k\n,0) = 

A{n - k,k\n,0) = 

A{k,n-k\n,0) = 

A{n - k,k\n,0) = 




2.2 Two Input Beams 

We now calculate the general probability P{Ni,ni + n2 — A''i|ni, 77-2) that A''! output photons 
are observed along the direction of input beam 1 when the number of photons in the input 
beams in ni and n2- 



We first give a classical wave analysis. The input waves have amplitudes 01^2 = yni^2fi^, 
and are in phase at the center of the beam splitter. The output amplitudes are the sums 
of the reflected and transmitted parts of the input amplitudes. A relected amplitude has a 
phase shift of 90° relative to its corresponding transmitted amplitude, as discussed in sec. 2.1. 
In the 50:50 beam splitter, the magnitude of both the reflected and transmitted amplitudes 
from a single input beam are times the magnitude of the amplitude of that beam. 

Hence, the output amplitudes are 

Ai = -^(01 + ^2), (12) 

A2 = -^(^1 + 02). (13) 

Taking the absolute square of eqs. p2|l - (|T!^ . we find the output beams to be described by 

iVi2 = ^ = ^^ = ^^^. (14) 
hu 2huj 2 
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The classical view is that a 50:50 beam splitter simply splits both input beams, when they 
are in phase. 

For a quantum analysis, we proceed by noting that of the A'^i photons in output beam 1, 
k of these could have come by transmission from input beam 1, and Ni — khj reflection from 
input beam 2 (so long a.s Ni — k < 77-2). The probability amplitude that k out of A'^i photons 
are transmitted from beam 1 while Ni — k photons are reflected from beam 2 is, to within 
a phase factor, the product of the amplitudes for each of these configurations resulting from 
a single input beam: 

; / 1 \ ("i+"2)/2 

A(k, N, - k\m, 0)A{N, -k,n2-N, + k\0, n,) = {-ir-'^JC^'C^^k [^) , (15) 

referring to eqs. (|H|l- |TT|) . The most dramatic nonclassical features to be found below can be 
attributed to the presence of the factor (— that arises from the 90° phase shift between 
reflected and transmitted photons. 

Since photons obey Bose statistics, we sum the sub-amplitudes (fTSj) . weighting each one 
by the square root of the number of ways that k out of the A'^i photons in the first output 
beam can be assigned to input beam 1, namely C^^, time the square root of the number 
of ways that the remaining ni — k photons from input beam 1 can be assigned to the 
photons in output beam 2, namely C^^_^ to obtain^ 

AiNi,ni + n2- Ni\ni,n2) = E \/Ck'C^'-kMk, iVi - k\ni, 0)A(iVi -k,n2-N^ + k\0, n^) 

k 




When evaluating this expression, any binomial coefficient in which m is negative, or 
greater than n, should be set to zero. 
The desired probability is, of course, 

P{Ni,ni + n2 - Ni\ni,n2) = \A{Ni,ni + n2 - Ni\ni,n2)\^ (17) 

Some examples of the probability distributions for small numbers of input photons are 
given below. 

2.2.1 Two Input Photons 
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Output (iVi, 


A^2| 


\ni,n2) 


(0,2| 
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|2,0) 
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1 
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1 
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^Delicate to justify not also including factors C^^_i^, and Cj^^_^^^_^j, these being the ways of assigning 
photons to output beam 2 - but these factors are the same as those already included, and so should not be 
counted twice... 
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When rii or 77,2 is zero, the probabihty distribution is binomial, as found in sec. 2.1. When 
ni = = 1 there is complete destructive interference between the cases where both photons 
are reflected (combined phase shift = 180°) and when both are transmitted (combined phase 
shift = 0). This quantum result is strikingly different from the classical expectation that 
there would be one photon in each output beam. 

2.2.2 Three Input Photons 
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Output {Ni,N2\ 
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2.2.3 Four Input Photons 
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2.2.4 Symmetric Input Beams: n\ = n2 = n 

There is zero probability of observing an odd number of photons in either output beam. 

To see this, note that when rii = 77,2 = n, the magnitudes of the subamplitudes are equal 
for having k photons appearing in output beam 1 from either input beam 1 or input beam 
2. However, the phases of these two subamplitudes are 180° apart, so that they cancel. In 
particular, when k photons are transmitted into output beam 1 from input beam 1, then 
Ni — k photons are reflected from input beam 2 into output beam 1; meanwhile, n—k photons 
are reflected from input beam 1 into ouptput beam 2. So the overall phase factor of this 
subamplitude is ^■'Vi-fc+^-fc = (—1)*' i""*"^^. Whereas, if k photons are reflected from input 
beam 2 into output beam 1, then Ni — k photons are transmitted from input beam 1 into 
output beam 1, and so n — Ni + k photons are reflected from input beam 1 into output beam 
2. So the overall phase factor of this subamplitude is = (—1)'^ The phase 
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factor between these two subamplitudes (whose magnitudes are equal) is 2^^^ = (— 1)^^, 
which is —1 for odd Ni, as claimed. 

For the case of observing an even number of photons in the output beams, a remarkable 
simplification of eq. (fTBj) holds jS] . I have not been able to show this by elementary means. 
It does follow by inspection when m = or n, in which case eq. I\W\) contains only a single 
nonzero term. In general, the index k in eq. ifl^) for A{2m,2n — 2m\n,n) runs from to 
2m if 2m < n, or from 2m — n to n if 2m > n. There are an odd number of terms, the 
central one having index k = m. By a strange miracle of combinatorics, the sum collapses 
to a simplified version of the central term of the series.... Namely, 

Ai2m,2n - 2m\n,n) = (-1)""™ (^) V^-™^--"- (18) 
Therefore, the n + 1 nonvanishing probabilities for symmetric input beams are 

. 1 s 2n -1 

P(2m, 2n - 2m|n, n) = - C^Cll-^^ ^ , , (19) 

y n ^ n ' 

where the approximation holds for large m and large n. Note that /q^ dxj \Jx{l — x) = vr. 
This probability distribution peaks for m = or ra, i.e., for all photons in one or the other 
output beam, with value 

P{0,2n\n,n) = P{2n,0\n,n) = (-] C^". (20) 

The probability of finding all output photons in a single beam when the input beams are 
symmetric is larger by a factor C^" than when there is only a single input beam (of the 
same total number of photons), because there are C^" ways of assigning the n photons from 
input beam 1 to the 2n photons in the output beam. This is an extreme example of photon 
bunching caused by the beam splitter. 

However, I reamin uncomfortable with the result iflPj) because it does not agree with the 
classical prediction l\14]) in the large n limit 



3 Appendix: Phase Shift in a Lossless Beam SpHtter 

We give a classical argument based on a Mach-Zender interferomter, shown in the figure 
below, that there is a 90° phase shift between the reflected and transmitted beams in a 
lossless beam splitter. Then, following Dirac's dictum, we suppose that this result applies 
to a single photon. 

A beam of light of unit amplitude is incident on the interferometer from the upper left. 
The reflected and transmitted amplitudes are re^'^'' and te*"^* , where magnitudes r and t are 
real numbers. The condition of a lossless beam splitter is that 

r^ + t^ = l. (21) 

The reflected and transmitted beams are reflected off mirrors and recombined in a second 
lossless beam splitter, identical to the first. 
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Then, the amphtude for transmission at the first beam sphtter, followed by reflection 
at the second, is tre^^^*~^'^^\ etc. Hence, the recombined beam that moves to the right has 
amplitude 



while the recombined beam that moves downwards has amplitude 

The intensity of the first output beam is 

= |Ai|'=4r¥, 

and that of the second output beam is 

/2 = 1^2!^ = r^ + t^ + 2rH'^ cos2(0f - 0^). 
For lossless splitters, the total output intensity must be unity, 

Ji + /2 = 1 = (r^ + ^2)2 + 2rV[l + cos2(0t - 0,)]. 
Recalling eq. (EH); we must have 

(Pt-(pr = ±90°, 

for any value of the splitting ratio r"^ -.t^. 



(22) 
(23) 

(24) 

(25) 

(26) 
(27) 
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